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Using a combination of first-principles simulations, based on the density functional theory, and 
Boltzmann's semiclassical theory, we have calculated the transport and thermoelectric properties of 
the half-metallic two dimensional electron gas confined in single SrRuOa layers of SrTiOs/SrRuOs 
periodic superlattices. Close to the Fermi energy we find that the semiconducting majority spin 
channel displays a very large in-plane component of the Seebeck tensor at room temperature, S = 
1500 /iV/K, and the minority spin channel shows good in-plane conductivity a — 2.5 (mQcm)~^. 
However, contrary to the expectation of Hicks and Dresselhaus model about enhanced global ther- 
moelectric properties due to the confinement of the metallic electrons, we find that the total power 
factor and thermoelectric figure of merit for reduced doping is too small for practical applications. 
The reason for this failure can be traced back on the electronic structure of the interfacial gas, 
which departs from the free electron behaviour on which the model was based. The evolution of the 
electronic structure, electrical conductivity, Seebeck coefficient, and power factor as a function of 
the chemical potential is explained by a simplified tight-binding model. We find that the electron 
gas in our system is composed by a pair of one dimensional electron gases orthogonal to each other. 
This refiects the fact the physical dimensionality of the electronic system can be even smaller than 
that of the spacial confinement of the carriers. 

PACS numbers: 72.20.Pa,73.40.-c,74.70.Pq,73.20.At 



I. INTRODUCTION 

The interest in thermoelectric materials has undergone 
a revival over the last decade.^ ^ The reason behind this 
resurgence is double. On the one hand, recent experimen- 
tal^ and theoretical^ discoveries allowed to glimpse 
new ways to significantly improve the efficiency of these 
materials, quantified by the thermoelectric adimensional 
figure of merit 



ZT = T, (1) 

where S is the Seebeck coefficient (also called ther- 
mopower and denoted by a by some authors), a is the 
electrical conductivity, Kq is the electronic (lattice) 
contribution to the thermal conductivity, and T is the ab- 
solute temperature. On the other hand, there is an indu- 
bitable technological interest over these systems, as sig- 
nificant improvements in this field will potentially impact 
fuel consumption efficiency and allow to build diminute 
cooling devices without moving parts. ^ 

It is usually accepted that applications in this field will 
only be cost effective for materials where ZT is signifi- 
cantly larger than 1. However, the task of increasing 
its value remains challenging, since all the participating 
material's parameters in Eq. ([T]) are strongly intercon- 



nected, and also dependent on material's crystal struc- 
ture, electronic structure and carrier concentration.!^ In- 
creasing the Seebeck coefficient for simple materials re- 
quires a decrease in the carrier concentration, that yields 
to a concomitant reduction in the electrical conductivity. 
Also, an increase in the electrical conductivity leads to a 
comparable increase in the electronic contribution to the 
thermal conductivity (as expressed in the Wiedemann- 
Franz law.) In conventional solids, a limit is rapidly ob- 
tained where a modification in any of these parameters 
adversely affects other transport coefficients so that the 
resulting ZT for a given material at a given temperature 
does not vary significantly.^^ 

Typical good thermoelectrics are doped semiconduc- 
tors with intermediate values of the carrier density (close 
to 10^^-10^^ cm~^) displaying, at the same time, large 
Seebeck coefficients and good charge mobility. A typical 
example^ of this behavior is Bi2Te3 that, after doping, 
displays a ZT value close to 1 at 320 K with a resitiv- 
ity ~ 1 ml^cm and a Seebeck coefficient of 225 /iV/K. 
Other heavy- metal-based materials, such as group IV 
chalcogenides (lead telluride, PbTe, and related mate- 
rials), exhibit large figure of merit at intermediate tem- 
peratures (up to 850 K). Mo re inte resting materials such 
as clathrates,^^ half-HeuslerLifl^ skutterudites,!^^!^ or 
strongly correlated oxides. ^^^^ have also been identified. 
Finally, i t is n oteworthy that graphene layers,Sl silicon 
nanowire J^^^ and molecular junctions have also received 
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much attention due to their low thermal conductivity. 

A completely new route to enhance the figure of 
merit was opened by Hicks and Dresselhaus in the early 
nineties/ With a theoretical model, based on a sim- 
plified electronic band structure (assuming free-electron 
parabolic bands in two dimensions and a localized bound 
state in the third direction), these authors showed that 
a significant increase in ZT would be possible due to 
the modification of the electronic properties of some ma- 
terials when prepared in the form of quantum-well su- 
perlattices ^ or nanowires.^ In particular they predicted 
that production of bidimensional heterostructures has 
two main consequences that improve the value of ZT\ 
(i) The confinement of the charge carriers in a plane is 
expected to reduce the dispersion of the density of states 
(DOS) of the bulk system, creating a sharper DOS at 
Fermi energy, and (ii) the nanostructuration of the sys- 
tem along a given spatial direction also favors the dis- 
persion of phonons that have wavelengths longer than 
the period of the superlattice, leading to a decrease on 
the thermal conductivity due to scattering of the lat- 
tice vibrations by the interface between layers. The first 
condition is in line with Mahan and Sofo's proposaP for 
an "ideal thermoelectric". From a purely mathematical 
point of view, these authors found how a (5-shaped trans- 
port distribution maximizes the thermoelectric proper- 
ties. 

Experimentally, these ideas have been thoroughly 
checked in semiconducting heterostructures containing 
heavy non-metal ions. As a proof of concept, pio- 
neering experimental works were carried out by Hicks 
and coworkers^^ in PbTe/Pbi_a^Eua^Te multiple quantum 
wells, showing a good agreement between the experimen- 
tal results and the theoretical model predictions for the 
increase in the figure of merit. Values of ZT as high as 
2.4 have been measured in p-doped Bi2Te3/Sb2Te3 su- 
perlatt ices .1^ However, these systems are far from ideal as 
they decompose at the temperature where thermoelectric 
materials are expected to function (they start to decom- 
pose at T ^ 200 C), and contain poisonous elements like 
lead or bismuth. 

Also, in good agreement with the Hicks and Dres- 
selhaus prediction, Ohta et found that periodic 
SrTiOs/Nb-doped SrTiOa superlattices, where a two- 
dimensional electron gas (2DEG) is formed at the Nb- 
doped layer, exhibit enhanced Seebeck coefficients as the 
width of the doping layer is reduced. In particular, when 
its thickness reaches the ultimate thickness of one unit 
cell a very high value of 5* = 850 /iV/K was observed, 
and a ZT of 2.4 for the 2DEG was estimated (this corre- 
sponds to an effectice ZT of 0.24 for the complete device 
having the 2DEG as the active part). However, this en- 
hancement of the thermoelectric properties in 2DEG at 
oxide superlattices seems not to be so universal. Recent 
experimental result^ could not find any enhacement of 
the Seebeck effect due to the electronic confinement in 
the metallic state at the n-type LaAlOa/SrTiOs inter- 
face.!^ 
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FIG. 1. (Color online) (a) Schematic representa- 

tion of the unit cell, periodically repeated in space, of 
(SrTi03)5/(SrRu03)i superlattices. Ti atoms (in blue) Ru 
ones (in gray) are situated at the center of the octahedra, as 
denoted by the layer labeling, while O (in red) is placed at 
the vertex and Sr (in green) at the interstices, (b) Layer- by- 
layer PDOS on the atoms at the SrBOs {B = Ti or Ru) for 
the corresponding layer at the same height as in (a). Major- 
ity (minority) spin is represented in the upper (lower) half of 
each panel. 



Very recently, we have proposed an alternative mech- 
anism to generate 2DEG at oxide interfaces, playing 
with the possibility of generating a quantum well based 
on the different electronegativity of the cations in the 
perovskite structures used to build the interface. In par- 
ticular a half-metallic spin-polarized 2DEG was theoret- 
ically predicted in (SrTi03)5/(SrRu03)i superlattices. 
The electron gas is fully localized in the SrRuOa layer 
due to higher electronegativity of Ru^+ ions compared 
to Ti'^^ ones. The 2DEG presents a magnetic moment of 
\i = 2/iB/Ru ion. All these properties can be seen in the 
layer-by-layer projected density of states (PDOS) shown 
in Fig. [l] where the only contribution to the DOS at the 
Fermi energy comes from the minority spin SrRuOa layer, 
while the majority spin channel behaves like a wide-gap 
semiconductor. From the point of view of the design of 
thermoelectric materials [see Eq. ([T])] this system seems 
very promising due to the possible combination of a high- 
Seebeck coefficient (coming from the semiconducting ma- 
jority spin channel) and metallic conductivity (coming 
from the minority spin channel). 

The main objective of this work is to assess the thermo- 
electric properties of the 2DEG in (SrTi03)5/(SrRu03)i 
superlattices, and to ascertain whether the Hicks and 
Dresselhaus model is applicable or not in this kind of 
systems. In order to achieve this goal, and starting from 
accurate first-principles electronic structure simulations 
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on this interface, we use the Boltzmann transport the- 
Qr-^ |29 | 30| the constant scattering time approxima- 

tion to obtain the Seebeck coefficient and other transport 
functions. 

The rest of the paper is organized as follows. Com- 
putational details are summarized in Sec. [llj The main 
results on the Seebeck coefficient, electrical conductiv- 
ity, and power factor for the superlattice are presented in 
Sec. |III[ Finally, in Sec. |IV| we discuss a simple model to 
understand the physical origin of the previous transport 
coefficients. 



II. COMPUTATIONAL DETAILS 

We estimated the electrical conductivity and the See- 
beck coefficient through the semiclassical Boltzmann the- 
ory within the constant relaxation time approximation, 
as implemented in the BoltzTraP code.^^ This imple- 
mentation relies on the Fourier expansion of the band- 
energies, provided by a first-principles electronic struc- 
ture code. Following our previous work, we have used 
siesta!^ to compute both the relaxed atomic and elec- 
tronic band structures of (SrTi03)5/(SrRu03)i super- 
lattice in the local density approximation (LDA) to the 
density functional theory. An extra Hubbard- term, 
following the rotationally invariant LDA-\-U scheme of 
Anisimov et al.^^^ is included to account for the strong 
electron correlations, with a /7eff of 4.0 eV applied only to 
the d orbitals of Ru, as in Ref. 28, In order to get smooth 
Fourier expansion of the one-electron eigenenergies and 
converged transport coefficients we proceed in a two step 
procedure: (i) first we relax the atomic structure and the 
one-particle density matrix with a sensible number of k- 
points (12 X 12 X 2 Monkhorst-Pack mesH^, and (ii) 
freezing-in the relaxed structure and density matrix, we 
perform a non-self-consistent band structure calculation 
with a much denser sampling of 72 x 72 x 17 (5994 /c-points 
in the irreducible Brilloin zone; 88128 /c-points in the full 
Brillouin zone). The rest of the computational param- 
eters remain the same as in Ref. 28, The robustness of 
the results presented below have been doubled-checked 
using the CRYSTAL09 cod^ within the Bl-WC hybrid 
functionaP^ that mixes the generalized gradient approx- 
imation of Wu and CoherP^with 16 % of exact exchange 
with the Bl scheme.^^ 

The atomic structure of the (SrTi03)5/(SrRu03)i su- 
perlattice at low temperature includes the rotation of 
both the TiOe and RuOe octahedra along the tetragonal 
z-axis of the system [to establish the notation, we will 
call the plane parallel to the interface the (x^y) plane, 
whereas the perpendicular direction will be referred to 
as the 2;-axis]. Taking into account that the tempera- 
ture at which bulk SrTiOs undergoes the tetragonal to 
cubic transition is only 105 K, we expect that these dis- 
tortions are fully suppressed when the system is acting 
as a thermoelectric generator at room or higher temper- 
atures. Therefore, in the present study we do not allow 



rotation and tiltings of the oxygen octahedra during the 
atomic relaxations of the superlattices. 

With the first-principles band structures, e^^k, com- 
puted as indicated above, together with the space group 
symmetry of the superlattice, we feed the BoltzTraP 
code. There, after performing the Fourier expansion, the 
conductivity tensor can be obtained as 



crc,/3(i,k) = eVi,k^a(^,k)v^(i,k), 



(2) 



where e is the electronic charge, r^^k is the relaxation 
time, and 



Idsi 



% dka 



(3) 



is the a component of the group velocity for an electron 
in band i. Now, from the previous conductivity matrix 
we can compute the relevant transport tensors that re- 
late the electric current with an external electric field 
[aa^{T,fi)] or temperature gradients [z^q,/3(T, /i)]. These 
tensors depend on the temperature, T, and the chemical 
potential, /i, that determines the number of carriers or 
the level of doping. The final expressions are given by 



r dk 




1 8^ 


de 



crc,/3(^,k), (4) 



and 



r dk 


\ df{T,ii^ 


1 8^ 





aai3{i,k)[s{k)-ij], 
(5) 

where / is the Fermi-Dirac distribution. Finally the com- 
ponents of the Seebeck tensor can be computed as 



(6) 



The electronic contribution to the figure of merit is 
summarized in the power factor, PF = S'^a, which is 
the numerator of the right-hand side in Eq. ([T]). For a 
magnetic system the value of the PF can be calculated 
from the individual spin bands using Eq. (|4| and Eq. ^ 
and 



PF 



While computing the previous transport properties two 
major approximations are considered: (i) the relaxation 
time r is treated as constant, independent of tempera- 
ture, band number, occupation and the k vector direc- 
tion. While this approximation is fairly strong, tests of 
this method^^ with semiconductor thermoelectrics like 
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Bi2Te3 and systems including electron correlation like 
CoSba lead to a reasonable agreement with experiment. 
After Eq. (|6|, this approximation allow us to compute 
the Seebeck coefficient on an absolute scale (independent 
of r). However the conductivity can be calculated only 
with respect the relaxation time, and a value of r has 
to be introduced as a parameter [typically using the the- 
oretical cr/r value obtained from Eq. Q to reproduce 
exactly the experimental conductivity at a given temper- 
ature and carrier density n]. Here, the relaxation time 
value r = 0.43 x 10~^^ s employed in the calculations was 
obtained from fitting the room temperature conductivity 
a = 1.667 X 10^ S/m of bu lk SrT iOs at electron concen- 
tration n = 1 X 10^^ cm~^,'^^^which is very similar to 
that obtained from SrRuOs .■^^EH 

The second approximation is the "rigid band ap- 
proach" that assumes that the band structure does not 
change with temperature or doping, and therefore is fixed 
independently on the chemical potential. 



III. RESULTS 
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FIG. 2. (Color online) Total DOS (solid black lines) and 
PDOS on some Ru(4d) orbitals, showing the main character of 
the bands around the Fermi energy in (SrTi03)5/(SrRu03)i 
superlattices. Green dot-dashed lines are the projections on 
the 4dxy orbitals, blue dashed lines on the 4:dxz,yz orbitals 
(these curves are exactly degenerate) , and red dotted lines on 
the 4:d^2_y2 orbitals. The zero of energies is aligned at the 
Fermi level. The wiggles around zero are caused by the finite 
k resolution. 



A. Electron localization and band structure 

In Fig. [l] we show the DOS of the 

(SrTi03)5/(SrRu03)i superlattices projected layer- 
by-layer. As previously discussed in Ref. |28l the 
electronic structure displays a half-metallic state where 
the conduction only takes place for the minority spin 
electrons that are completely confined to the SrRuOa 
layer. For the majority spin our LDA+/7 calculations 
predict a gap of ^ 1.4 eV around the Fermi energy, 
which is typical for wide-gap semiconductors. 

A more detailed description of the nature of the states 
around the Fermi energy is obtained when we plot the 
PDOS for the M orbitals of Ru (Fig. ^. Since Ru'^+ is a 
low-spin transition metal ion, we expect the conduc- 
tion band to have a strong t2g{4:dxyAdxzAdyz) charac- 
ter. Indeed, the conduction band in the minority spin is 
formed by the half- filled degenerate Ru{4:dxz,yz) orbitals 
while the Ku{4:dxy) band is mainly situated at an energy 
slightly above the Fermi level. Similarly, the majority 
spin valence band is composed by Ru t2g orbitals, but 
since they are shifted to lower energies they are strongly 
mixed with 0{2p) bands. In Ref. 28 these features were 
explained with the use of a simplified tight-binding model 
including three main physical ingredients: (i) the bidi- 
mensionality of the SrRuOs layer, (ii) the difference of 
in-plane and out-of-plane bonding for the Ru-ions, and 
(iii) the electron-electron interactions as described by a 
Hubbard term. In Fig. [2] we can also see the Ru(4(ia,2_^2) 
band, however as it lies at relatively high-energies its 
presence is negligible when discussing the transport prop- 
erties in this system. 



B. Transport calculations 

In Fig. |3] we show the calculated electrical conductiv- 
ity, Seebeck coefficient and power factor for T = 300 K as 
a function of the position of the chemical potential (i.e. 
doping level) for both the majority and minority spin 
channels. We also compare them to their corresponding 
DOS. Due to the tetragonality of the superlattice, all the 
previous transport tensors are diagonal with only two in- 
dependent components: one parallel to the interface {xx 
= yy)^ and a second one perpendicular to the interface 
{zz component). Since all the carriers are confined to 
move in the SrRuOa plane, from now on we will only 
focus on the parallel one. 

As expected for a half-metal, the behavior of these 
quantities around the Fermi energy at zero doping is very 
different according to the different nature of each spin 
channel. 

In the case of the minority spin, the system is metal- 
lic in the SrRuOa layer (Fig. |2|, and the conductivity 
presents a local maximum close to /i = 0. Then, a de- 
creases quickly as the chemical potential decreases and 
gets closer to the Ru {4:dxz,yz) band edge (around -0.6 
eV below the Fermi energy, see Fig.|2]). The other spm 
component (the majority one) is semiconducting and the 
small non-zero contribution to a (indiscernible from zero 
in the scale of the figure) comes from the use of a finite 
temperature in the simulation (electrons thermally ex- 
cited to the conduction bands, leaving behind holes in the 
valence band) . In both cases the conductivity quickly in- 
creases when > 1.5 eV as the doping starts to involve 
the large density of states associated with SrTiOs levels. 



5 




-2-1012-2-1012 




Chemical potential (eV) Chemical potential (eV) 



FIG. 3. Results of the calculation of the electronic DOS and 
transport properties at 300 K [conductivity (cr), Seebeck co- 
efficient (S), and power factor (PF)] of (SrTi03)5/(SrRu03)i 
superlattices as obtained by Boltzmann's semiclassical trans- 
port theory. On the left column we show results for the ma- 
jority spin bands and on the right for the minority spin ones. 



When we compare these values with those of typical high- 
efficiency thermoelectrics, a ~ l{mQcm)~^ ^ we find that 
the value for the minority spin is larger 3{mQcm)~^] 
while that of the majority channel is much smaller. 

Similarly, the curves for S also refiects the different 
nature of the spin components. In the case of the major- 
ity spin we note a discontinuity of the Seebeck coefficient, 
showing the change from hole (region with /i < 0) to elec- 
tron (region with /i > 0) doping. The maximum absolute 
value for this channel is 5* ~ 1500 /iV/K and occurs for 
relatively small hole dopings {/a = —0.25 eV) This value 
is significantly larger than the room temperature S — 
480 /iV/K found by Ohta et al? for the 2DEG in one 
unit-cell-thick Nb-doped SrTiOa superlattices, than the 
value of 5 = 1050 /iV/K in TiOs/SrTiOs heterointerfaces 
with an electron concentration of 7 x 10^^ cm~^,'^ and 
also than those typical associated to good bulk thermo- 
electrics, which are usuall^^ around 150-250 /iV/K (the 
reader has to keep in mind that we are comparing maxi- 



mum values of the Seebeck coefficient, thought they can 
be achieved at different carrier concentrations of holes or 
electrons.) The opposite happens to the minority spin 
channel, where the Seebeck coefficient S is very small at 
/i = as the conduction changes from being dominated 
by electrons to holes in the semi filled Vc^{^dxz,yz) bands. 
When the system is doped, the absolute value increases 
linearly with the chemical potential but the energy scale 
is two orders of magnitude smaller than in the majority 
spin case, a typical factor when comparing the Seebeck 
coefficients of metals and semiconductors. Only when 
the chemical potential is close to the lower edge of the 
conduction band (/i ~ — 0.6 eV) the minority spin shows 
a pronounced enhancement of the thermopower reaching 
a moderate value of 5 = -35 \i V/K. 

Calculation of the power factor for each of the spins 
shows strong compensation of Seebeck coefficient and 
conductivity in both channels giving rise to very small 
values (< 0.5 /iW/K^cm) around the Fermi energy. In- 
deed, while the majority spin displays a large Seebeck 
coefficient and a very small conductivity, in the minor- 
ity spin a reversed situation is found. In both cases the 
power factor for the undoped system is almost negligible 
when compared to those of good thermoelectrics (20 - 
50 /iW/K^cm). Only for very strong hole doping, close 
to one hole per Ru^+ ion, an appreciable enhanment is 
observed for the power factor, both for the majority spin 
component, PF=3.6 /iW/K^cm and the minority spin 
one, PF=1.7 /iW/K^cm. This result is summarized in 
Fig. |4| where the total PF is calculated as indicated in 
Eq. Q. There we can observe that the total power fac- 
tor is still very small for low dopings, and only when the 
system is strongly hole-doped (>le/Ru) an appreciable 
PF is obtained. 
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FIG. 4. Total PF for the (SrTi03)5/(SrRu03)i superlattices 
at 300 K as a function of the chemical potential. 
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IV. DISCUSSION 

In the previous section we have found that even though 
(SrTi03)5/(SrRu03)i superlattices exhibit a 2DEG, 
whose width is confined just within a single perovskite 
unit cell, the power factor does not display the strong en- 
hancement expected by Hicks and Dresselhaus7 In order 
to understand the previous results we will analyze them 
using both a free electron and a tight-binding model. The 
choice of the free electron model was motivated by the 
fact that the Hicks and Dresselhaus works were based on 
parabolic bands to describe the electronic structure of 
the system, while the tight-binding model has proven to 
be very successful to explain the first-principles results of 
the superlattice under study. 

We define a simple one band free electron model in one, 
two or three dimensions using the energy dispersion, 

i 

where is the (isotropic) effective mass and i runs over 
the dimensions of the system. Similarly, we define the 
energy dispersion for a one-band tight-binding model as 

e{k) = 27^cos(A:ia), (9) 

i 

where 7 is the characteristic interaction energy (band 
width) of one of the tight-binding center with its first- 
neighbors along (100), and a is the cubic lattice spacing. 
A graphical representation of the free electron and tight- 
binding bands is presented in Fig. |5] 
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FIG. 5. (Color online) Graphical representation of (a) one 
dimensional tight-binding, and (b) free-electron bands that 
illustrate, respectively, Eq. ^ and Eq. ([8|. In (a) we show 
with red dashed parabolae the free-electron band approxima- 
tion to the tight-binding model at the band-edge while the 
blue dash-dotted line depicts the group velocity at mid-band 
{e = 0) and the band-edge {e = ±27). 

In Fig. [6] we compare the DOS and transport 
properties for both models using the formalism de- 
veloped in Sec. [Hj For ease of comparison with the 
(SrTi03)5/(SrRu03)i superlattice, this plot is struc- 
tured in a similar way to Fig. [Sj 



The DOS plots correspond with the textbook exam- 
ples'^ of the corresponding models where, for exam- 
ple, the 2D free-electron model involves a step function. 
Comparing both models for the same dimensionality we 
observe that, at the band edge, both are quite similar 
with a logarithmic divergence (in ID), a finite discon- 
tinuity (in 2D), and a functional dependence with the 
energy cx ^/e (in 3D). However, the models quickly differ 
at higher energies. In particular, the tight-binding bands 
are symmetrical around the center of the band, while the 
free-electron ones are not. The reason behind the simili- 
tudes between the two models comes from the fact that 
the tight-binding dispersion curves can be approximated 
by parabolae (see Fig. [5| at the band-edges. Compar- 
ing these simple models to the first-principles calculated 
DOS of the conduction band (top-right panel of Fig. [3| 
we find that the most similar one is the ID tight-binding 
model. In particular, both pictures display a very char- 
acteristic two-peaked structure. While it is reasonable 
that a tight-binding model is more adequate than a free- 
electron one to describe the narrow 4(i-bands of Ru, it 
seems surpring that the DOS resembles that of a IDEG 
rather than a 2DEG one. The reason for this behavior 
is that in the case of Ru{4:dxz) and Ru{4:dyz) bands the 
hopping parameter in the conducting plane is only large 
along X or y directions, respectively [see Fig.[7|;a)]. Thus, 
the 2DEG in the (SrTi03)5/(SrRu03)i superlattice is, in 
fact, formed by two half- filled orthogonal ID bands. In 
contrast, in the Bj\i{Adxy) and Bj\i{Adx'2-y'2) bands the or- 
bitals bond equally well in the x and y directions [see 
Fig.[7|^b) and[7]^c)] forming a proper 2DEG. Indeed, if we 
compare the DOS of the majority spin BM{Adx2_y2) band 
in Fig. [2] we observe a very similar shape to that of the 
ideal bidimensional tight-binding model shown in Fig. [6] 
For the BM{Adxy) a good agreement with the model can 
be also be achieved if the tight-binding expansion is ex- 
tended to include interactions with in-plane neighbors 
along (110) directions that shift the central DOS peak to 
higher energies. 

Regarding the electric conductivity, the results ob- 
tained for both the free electron and the tight-binding 
models are equivalent at the band-edge where the tight- 
binding bands can be approximated by parabolae. This 
can be seen in the way the conductivity curves decay in a 
quicker way as the energy gets closer to the lower bound 
of the band. However, the behavior for half filling is quite 
different in both models. In particular, the tight-binding 
model predicts a maximum conductivity in the middle 
of the band, a behavior not observed for the free elec- 
tron approximation. This conductivity maximum corre- 
sponds to the maximum of the group velocity at £ = 
as deduced from Eq. ([9| and can be graphically visual- 
ized as the slope of the band diagram in Fig. [5] This 
maximum can clearly be seen around /i = in the full 
conductivity calculation for the minority spin channel in 
(SrTi03)5/(SrRu03)i (Fig.[3|, and the ID character of 
the band can also be observed in the abrupt reduction of 
conductivity around /i ~ —0.5 eV. 
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FIG. 6. (Color online) Density of states (DOS), electrical 
conductivity (a), Seebeck coefficient (changed sign, —S) and 
power factor (PF) for a free electron model as described in 
Eq. (|8} (left panels), and the tight-binding model of Eq. ^ 
(right panels). Different colors represent the dimensionality 
of the system: solid black lines for three-dimensional, blue 
dashed lines for two-dimensional, and red dash-dotted lines 
for one-dimensional. 



The Seebeck coefficient follows a similar trend, with 
the free electron and tight-binding results closely match- 
ing each other at the edge of the band. At this point, the 
absolute value of the Seebeck coefficient is larger for bulk 
(3D) than in lower dimensionality systems (2D, ID). This 
may come as a surprise when it was experimentally shown 
that in the 2DEG in Nb-doped SrTiOj^the Seebeck con- 
stant increases with reduced dimensionality. However in 
that case it was argued^ that polarons are responsi- 
ble for transporlP3] and thus the present model cannot 
account for their behavior. Going back to the compari- 
son between both the free electron and the tight binding 
models we find that in the former case electronic bands 
are either for electrons or holes depending on the sign of 
the curvature in Eq. ^ around A: = 0, while in the latter 
the bands transit from electron {S < 0) to hole {S > 0) 
dominated with filling. Thus, the value 5* = is achieved 



y" 

► 

X 

FIG. 7. (Color online) Illustration of the Ad orbitals of Ru 
in the Ru02 layer and hopping parameters, 7, for (a) the dxz 
(an equivalent schema, rotated 90°, could be made for the 
dyz), (b) the dxy, and (c) the dx2_y2 orbitals. In (a), the 
hopping along z is strongly hindered as it involves moving 
into Ti levels which are much higher in energy. From the 
overlap of the orbitals it is clear than 7 >> 7^ Meaning of 
the colours for the atoms as in Fig. Ilja). 



exactly in the middle of the band, as observed both for 
the model in Fig. |6] and the full calculation for the mi- 
nority spin bands in the full heterostructure (Fig. |3|. 

Finally, we analyze the power factor in both mod- 
els. We observe that at the band edge the PF reaches a 
maximum value that is significantly enhanced when the 
dimensionality is reduced, in agreement with parabolic 
band models within effective mass theory used by Hicks 
and Dresselhaus.^^ Taking into account that PF=6'^<j it 
is surprising to find this enhancement knowing that the 
maximum value of the Seebeck coefficient is reduced with 
dimensionality. However, this decrease is over compen- 
sated by the increase of the conductivity at the edge of 
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the band. As the decay of the one dimension conductiv- 
ity is sharper close to the beginning of the band for lower 
dimensionality systems, the value of the conductivity is 
also larger close to the edge than in 3D leading to a global 
enhancement of the PF. In the middle of the band the 
tight-binding model predicts S = which leads to a null 
value of the PF. Thus, the fact that the conducting band 
is exactly half filled in (SrTi03)5/(SrRu03)i superlat- 
tices explains their very poor thermoelectric properties 
even though they display a half-metallic 2DEG with a 
width of a single SrRuOa layer and a very large ther- 
mopower. 

Therefore, the strong confinement of the electron 
gas in (SrTi03)5/(SrRu03)i superlattices does not in- 
volve a corresponding enhancement of the thermoelec- 
tric properties. Recent experiments in the conducting 
LaAlOa/SrTiOs interface also point to the absence of an 
enhancement of the Seebeck coefficient due to electronic 
confinement.^^ This seems to be at odds with the Hicks 
and Dresselhauss's model. ^ ^ However, there is no such 
discrepancy: simply, the physical approaches on which 
the previous mode was based are not fullfilled in the 
present case. The main difference of our system with 
those of Hicks and Dresselhauss is the fact that the bands 
in our system are narrow and half-filled, mostly due to 
the potential wells arising predominantly from the ionic 
charges, breaking the parabolic band approximation 
used in their calculations. 



V. CONCLUSIONS 

In this work we have theoretically studied the trans- 
port properties of the 2DEG present in the half-metallic 
(SrTi03)5/(SrRu03)i superlattice. Using Boltzmann's 
transport theory we have shown that: (i) the semicon- 
ducting spin channel displays a large Seebeck coefficient, 
S ~ 1500 /iV/K, larger than that found in Nb-doped 



SrTiOs thin films, (ii) however the total power factor is 
too small for the system to be a good thermoelectric ma- 
terial. The reason behind the low power factor is the half 
filling of the metallic spin channel that, while providing 
a high conductivity, finally quenches the total Seebeck 
constant of the system. 

Our results do not contradict those of Hicks et al^^ 
who predicted an important enhancement of the thermo- 
electric properties for systems whose transport is strongly 
confined in one or two dimensions although it shows its 
limits of application. Hicks and Dresselhaus assumed a 
semiconducting superlattice with parabolic bands that 
could be doped to increase the carrier density. How- 
ever, in our system the main transport bands are nar- 
row as corresponds to those with a strong Ilu{4:dxz,yz) 
character and are half filled being, as a consequence, not 
suitable for the application of the free electron model. 
Indeed, they are well described within the tight-binding 
approximation. Using this model we show that the 2DEG 
in (SrTi03)5/(SrRu03)i is unusual in the sense that is 
composed of two orthogonal bands where the hopping 
parameter is only strong in one dimension having thus, 
the properties of a IDEG. 
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